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The uniform electron gas, the traditional starting point for density-based many- 
body theories of inhomogeneous systems, is inappropriate near electronic edges. 
In its place we put forward the appropriate concept of the edge electron gas. 
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Since the work of Thomas and Fermi (TF) through 
the earliest papers on density functional theory (DFT) ||] 
up to ongoing work involving density gradients j^, the 
uniform electron gas has been the starting point of 
density-based approximate theories of inhomogeneous 
systems. However, real physical systems have electronic 
edge regions where effective single particle wave functions 
evolve from oscillatory to evanescent, and where clearly 
approximations starting from a uniform electron gas are 
ill-founded. 

The edge surface S, for a given system, may be pre- 
cisely defined by 



(1) 



where Veff{r) is the exact self-consistent potential of the 
Kohn-Sham (KS) equations and /i is the chemical po- 
tential, half-way between the highest occupied and the 
lowest unoccupied KS orbital energies (In what follows 
we take fi — 0). Outside of S all occupied KS orbitals 
decay exponentially (See Fig. 0). 

Much interesting science involves precisely such edge 
regions (See Fig. |^), such as ionization energies of the 
outermost electrons, molecular binding, surface energies 
etc. For these regions, in place of the uniform electron 
gas, we discuss in this paper the new concept of the edge 
electron gas, which is valid in the edge region near S . 

The edge gas concept is based on the principle of 
"nearsightedness" put forward in : The local electronic 
structure near a point r, while strictly speaking requiring 
a knowledge of the density n{r'), or equivalently of the 
effective potential Veff{r'), everywhere, in fact is largely 
determined by Veff{r') for r' near r. Nearsightedness 
applies just as much to the edge region as to the interior 
region far inside S. 

This principle has been applied to the density n{r) and 
to the one particle density matrix j^]. In the present pa- 
per we also apply it to the exchange energy as follows: 
We write 



n{r)R^ ^{r) dr . 



(2) 



where ^ ('''), the inverse radius of the exchange hole, is 
defined by 



R-\r) . - / !^dr' 



and Uxir; r') is the exchange hole density which satisfies 
/ n^{r;r')dr' = -1. 
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FIG. 1. Interior and Edge regions of a bounded system. S 
is the nominal dividing surface ( ). S { ) is the phys- 
ical dividing surface with I = jl. I is defind in Eq. ^ and 
1< 7 < 3. 

For a consideration of the electronic structure at a 
point r near S we first drop a perpendicular to the near- 
est point Te on S, take as origin of coordinates and the 
2;-axis through r (See Fig. |l|.); thus r = (0, 0, z). Near 
we have, up to first order in {x, y, z), 



Veff{r)~—Fz. (Airy Gas) 



(4) 



(3) 



This leads to the concept of the Airy gas (AG) , the sim- 
plest version of the edge gas (Fig. ^). The AG is a 
principal subject of this paper. Of course only the prop- 
erties of the AG near the edge at z = are of physical 
interest. 



FIG. 2. a) The Airy gas; a strictly linear v^ff, with a hard 
wall at the distant point L. b) The edge gas: a typical Veff{z) 
near the edge at 2 = 0. All levels up to /i = are occupied. 

In many cases Ve// varies rapidly along z but slowly 
along X and y, i. e. kI ~ 1, where I is the nominal surface 
"thickness" , defined in Eq. (^ below, and k is the curva- 
ture of S at Te- This leads to the more general concept 
of the Edge Gas (EG), 
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^eff{r)=veff{z); <^/(0) = -F, (Edge Gas) (5) 

(Fig. |b). The EG can be refined by including curvature 
corrections 0. 

There have been extensive previous studies of sev- 
eral model systems, notably by V. Sahni, A. G. Eguiliz 
and their collaborators 0, of electrons moving in effec- 
tive potentials Veff{z) with an electronic edge, say at 
z — 0. Their common features are: v^ff {~oo) > /i and 
?;e//(-t-oo) = constant < fi. The differences between ear- 
lier work and the present contribution are: (i) Our em- 
phasis on the locality of the electronic edge structure, i. 
e. its independence of the behavior of Veffif) far from 
the edge; (ii) our detailed analysis of the simplest uni- 
versal edge model, Veff{z) — —z (— oo < z < +00) 
which, by appropriate scaling, provides a first approxi- 
mation for dealing with general 31? systems with edges; 
(iii) preliminary discussions of further improvments. 

The Edge Electron Gas: We now consider the situation 
in Fig. ||b where the effective potential is independet of 
X and y but generic in z. 

The surface "thickness" / and the corresponding energy 
e are defined uniquely, up to constant factors, by 



\2mF J 



1/3 



2m 



1/3 



(6) 



where F = |</j:(0)|. 

The normalized eigenfunctions of the KS equations |^] 
are 

V'.(x,2/,z) = ^,(z)^e^('=^-+'=^^); i^„^,,)=S,,, (7) 

where p = (j, /cijfe); kiLi = 2'Kmi {i = 1,2) and A = 
L1L2 is the cross-sectional area. The functions ^Pj{z) 
satisfy 

(-^d^+ ^-^/Z^^) - ^ (8) 



and the boundary conditions 

(pj{~(X)) = (y5j(+Oo) = 



(9) 



The ej are negative and, for each ej, (fci,^) satisfy the 
inequality 



The density of the EG is given by 



(10) 



(11) 



where the factor 2 accounts for the spin and where Wj is 
the weight factor accounting for the number of occupied 
wave- vectors (fci, fc2). 



The exchange hole density, nx{r;r'), is defined by 
l|pi(r;r')P 

nx[r;r ) 



2 n{r) 



(12) 
(13) 



e„<0 



Inserting Eq. (|^) in Eq. (^3|) and integrating over ki 
and k2 gives 



pi(r;r') = - V (^^ (2)93^(2') 

TT ^ — ' 



kjJi{kjp') 



P' 



(14) 



where Ji (x) is the first order Bessel function and where, 
without loss of generality, we have set r = (0, 0, z), and 
have expressed r' in cylindrical coordinates p', a' and z' . 
kj is the maximum transverse wave number associated 
with ipj{z) given by kj = {2m\ej\/h'^)^^'^ . This leads to 

j j' 

x{Azy^ g{kjAz,kj,Az) , (15) 

where Az — \z — z'\ and the universal function g{s, s'), 
independent of Veff{z), is given by 

g{s,s') = ss' / dt 'I Lll^ ' . (16) 
Jo tvl + t 

We have calculated ^(s, s') numerically and present a 
contour plot in Fig. |^. Further details about g(s^ s') are 
available from A. E. Mattsson. 
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FIG. 3. The universal function g{s,s') of Eq. (|l|). 

The Airy Gas: We now consider the simplest model of 
an edge electron gas, the Airy Gas (AG) of Fig. ||a, in 
which the potential has the strictly linear form 



2 



. , J ~Fz ~oo < z< L {F>0) 
The normalized eigenfunctions ifjiz) satisfy 



(17) 



(18) 
(19) 



In suitable variables all the fj{z) are proportional to the 
Airy function Ai{—z); Ai(z) satisfies 



{-— + z)At{z) = 



Ai{+oo) = . (20) 
The function ^(z) = Ai{—z) is shown in Fig. ^. 




FIG. 4. The function $(2) = Ai{~z) 
Relevant properties of $(2) are S: 



$'(0) 
1 



$(0) = 0.35503; 
z ^ -00 : $(z) w ^^i/2|^|i/4 e' 
z ^ : <i>(z) « ^T7W sin (|NP/' + f ) 



0.25882 

' 2 1 



(21) 



In the AG all 3D states with energy up to /i = are 
occupied (Fig. ^). The AG is completely characterized 
by the length I and the energy e in Eq. (^). In units of I 
and e all AG's are identical. We shall later verify that for 
z » I the AG behaves like a locally homogeneous gas 
whose density is n{z) = {'iTi'^)-'^{2mF /fi^ f/'^ z'^/'^ . Thus 
the "thickness" of the edge region is ~ I. The fj{z) 
contributing to the edge electron density have Cj of the 
order — e, much less, in absolute value, than the lowest 
eigenvalue (^ —FL). These small e^, obtainable from 
the boundary condition fj{L) = 0, and the ^j{z) are 
given by 



-4^ 



1/2, 



ipj{z) 



rl/2 



$(- 

{Ll)~^/^ ^l 



(22) 



We note that the ej are equally spaced and that suc- 
cessive (^j(z)'s are identical except for a small shift, 
1{1/LY/'^tt, to the right. 

Using Eqs. (^) in Eq. (fl]) we obtain for the density 



n{z) 
noiC) 



(23) 
(24) 
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FIG. 5. The density of the Airy gas, exact ( ) and in the 

Thomas Fermi approximation ( ). 

Fig. H shows the exponential decay for C, « —1. For 

C»i 



no 



(0 = Ac^/^-T^r^/'sin(V2) + ... . (25) 



16n 



The first term is identical with the Thomas Fermi result 
(See Fig. ||); the second describes Friedel oscillations, due 
to the edge at z = 0. 

Exchange: Inserting the appropriate expressions for 
the AG, Eqs. (||), into Eq. (|^) gives 



R-\z)^r'R-o\0 , 



(26) 



where 



TT no(C) 

00 00 

J dC jdej de'$(C - e)$(C' - e)$(C - e')^{C " e') x 

xic-cr'5(%/^ic-c'i,v7ic-c'i) . (27) 







In Fig. H we see the exact R^^iO for the AG compared 
with its local density approximation (LDA) result 



(^x-o(C)) 



LDA 



3/2 (3A)V3 {no{Q)f' 



(28) 
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FIG. 6. The inverse radius of the exchange hole of the Airy 
Gas: ( ) exact, ( ) LDA. 

As expected the LDA exchange is much too small for 
z/l < —1, but accurate for z/l > 2. For ( — oo, 

For the AG the exchange hole density has been calcu- 
lated numerically and is shown for r = (0, 0, 0) in Fig. (|^). 
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FIG. 7. The exchange hole density nx{r') of the Airy Gas 
for r = (0, 0, 0) (see the black dot). The vertical line indicates 
the edge at 2 = where Veff[z) = ^l. The values on the 
contours are given in units of 10~^ Z^^. 

This may be compared with the isotropic exchange 
hole density of a uniform electron gas of density n, cen- 
tered at the point r and given by 



n^(r; r') = — 9n (sinr — r cost)^/(2 r®) , 



(29) 



where r = (37r^ri)"'^^^|r — r'|. 

Application to the Edge Problem of Thomas Fermi 
Theory: For simplicity we first consider electrons with 
chemical potential /i in a spherical effective potential 
Veffir). The edge surface S" is a sphere of radius 
given by Veff{re) = A*. The "edge problem" in TF the- 
ory is that~incorrectly-nTj;^(r) = for r > (Compare 
Fig. I). 



Following present ideas we consider points r = (0, 0, r), 
where r ~ re, and their reference point on 5, r^, = 
(0,0, Te). Near we approximate 

Veffir) ^^l-Fir,-r), F=Kff{r,)\ (30) 

which, by spherical symmetry, gives 

n{r) — l^^no{{re — r)/l) (edge region) (31) 

in the edge region which we take to be r > rg — with I 
defined in Eq. (^) and 7 = 2 (see Fig. |^) . In the interior 
region, r < re — 7/, we use the standard TF result 



n(r) 



3^(^(^ 



2m , 



3/2 



. (interior region) 
(32) 



7^ the two forms approximately agree (See 



At r = 
Fig. I). 

For an anisotropic Ue//(r), F{rf.) (= \yv{r)r=rA) and 
thereby I will vary over the surface S. 

Concluding remarks: The concept of "nearsightness" 
allows separate consideration of the interior and the edge 
regions. For the latter this paper begins to develop the 
concept of the edge gas, so far without correlation ef- 
fects. The joining of interior and edge regions, using the 
inverse of the exchange-correlation hole radius, R~}{r) 
(cp. Eq. 1^), will be described in a forthcoming work. 
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